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Abstract 

A  simple  solution  of  a  general  problem  involving  a  bifurcated 
•vravegtiide  is  presented.  The  purpose  of  the  work  is  to  explain  a  new 
and  simple  method  of  solving  such  problems  and  to  exhibit  an  organic 
connection  between  Sturm-Llouville  theory  aiod  the  theory  of  two-^)art 
boundary -value  problems. 
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1,         Introduction 

The  mathematical  aspects  of  the  theorj'  of  wave  propagation  in  longitudinally 
uniform  waveguides  are  discussed  in  Sturm-Li ouville  theory,  which  deals  with  the 
existence  of  a  set  of  experimentally  determinable  normal  modes  or  eigenfunctions 
and  corresponding  eigenvalues.     Regardless  of  the  transverse  variation  of  the 
electrical  properties  of  the   guide  in  particular  cases,  the  theory  furnishes  a 
list  of  qualitative  properties  which  the  eigenfunctions  and  eigenvalues  share 
with  all  other  eigenfunctions  and  eigenvalues  corresponding  to  the  same  boundary 
conditions,     ^rlhen  a  semi- infinite  bifurcation  is  introduced  into  the  guide,   two 
or  more   semi-infinite  waveguides  result;  the  difference  between  the  electrical 
properties  of  these  waveguides  is  mathematically  exhibited  in  a  change  of  boundary 
condition  or  interval  of  definition  of  the  modes  and  also  in  ths  new  eigenfunctions 
and  eigenvalues  that  ai-ise.     These  form  a  complete   set  of  functions  in  the  narrower 
waveguide  created  by  the  bifurcation.     The  qualitative  properties  of  this  new 
set  of  functions  are  not  the   same  as  those  possessed  by  the  functions  relating 
to  the  undisturbed  part  of  the  waveguide. 

The   present  paper  reports  an  investigation,    suggested  by  the  situation  des- 
cribed above,   of  the  possibility  of  expressing  the  solution  of  problems  relating 
to   the  bifurcated  waveguide  in  terms  of  the   solutions  relating  to  the  simpler  com- 
ponent waveguides.     Our  results  may  be  briefly  described  as  follows: 

It  is  found  possible  to  express  the  solution  of  the  two-part  problem  in  a 
very  simple  manner  in  terms  of  the  eigenfunctions  and  eigenvalues  which  relate  to 
the  component  waveguides.     Many  of  the  usual  complications  of  procedure  are  side- 
stepped,  and  the   solution  obtained  is  both  analytically  perspicuous  and  completely 
general.     The  proof  depends  purely  on  general  properties  of  Stum-Liouville  functions, 
i.e.,  on  the   qualitative  properties  of  the  modes  which  are   suitable  in  different 
portions  of  the  guide.     It  is  hoped  that  the  result  may  obviate   the  necessity  of 
subsequent  investigation  of  particular  cases,   and  that  it  may  make  it  possible  to 
assess  the  performance  of  a  bifurcated  variable-layer  waveguide  in  terms  of  the 
performance  of  unbifurcated  waveguides. 
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The  work  may  be  regarded  as  an  extension  of  an  earlier  Research  Report  L  J   in 
which  the  notion  of  two-part  boundary  value  problems  was  discussed.     However, 
familiarity  with  that  work  is  not  a  requisite   for  following  our  present  argument, 
whose   general  structure  is  as  follows. 

In  Section  2  we   formulate   a  waveguide  problem  dealing  with  a  symmetrically 
bifiircated  waveguide   supporting  one  incident  propagating  mode.     Problems  of  this 
kind  have   already  been  treated  by  various  authors  ••  ■'   for  the  case  of  a  homogeneous 
medium.     In  order  that  the  exposition  of  our  basic  idea  may  be  sufficiently  general, 
we   consider  the  case  of  a  medium  of  arbitrarily  varying  dielectric  constant.     How- 
ever, for  the  sake  of  simplicity,  we  assume  that  the  magnetic  permeability  is  con- 
stant,  and  that  the  dielectric  constant  is  symmetrically  distributed  with  respect 
to  the  bifurcation. 

In  Section  3  we  discuss  some  general  properties  of  the  solution  which  are 
guaranteed  a  priori  by  Sturm-Liouville  theory.     In  Section  U  we  construct  the 
solution  of  our  problem  on  the  basis  of  the  basically  well-known  information 
collected  in  Section  3.     Section  5  is  devoted  to  some  brief  remarks  on  the   re- 
lation to  Wiener-Hopf  theory. 

2.       Physical  origin  of  the  problem 

In  order  to  provide  a  motivation  for  the  mathematical  problem  to  be  consider- 
ed, we  call  attention  to  Fig.  1.     It  is  a  cross -sectional  representation  of  a 
symmetrically  bifurcated  parallel-plate  waveguide.     The  heavy  lines  represent 
the  conductors.     The  space  between  these  is  filled  with  a  medium  whose  magnetic 
permeability  \i  is  uniform  and  whose  dielectric   ■  constant"    ©(y)  is  even  in  the 
variable  y-a ,     In  this  medium  we  assume  an  electromagnetic  field  which  is  pro- 
pagated in  the  negative  x-direction,  with  a  frequency  co.     This  field  is  to  be  in- 
dependent of  2,   and  its  electric  vector  has  only  a  z-component.     Here  the  z- 
direction  is  taken  to  be  perpendiciilar  to  the  plane  of  Fig.  1.     We  assume  ft 
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time  dependency  of  the  form  e       • 
Then  from  Maxwell's  equations 
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Pfere  E  (x,y)  is  the   z-coraponent  of  the  field.     It  must  vanish  on  the  conductors, 
z 

Consequently,  if  we  denote  E    by  u,   and  use  the  fact  that  |j,  is  constant,  we  find 


that  we  must  solve  the  wave  equation 

2  2 

9  u  3  u     ^     ,  2,    . 

— ^     +     — ^     +     k   (y)u 

8x  9y 


0  <  y  <  2a, 


-00   <    X   <     00, 


with  boundary  conditions 

(2a)          u(x,0)   -  0, 

(2b)          u(x,2a)  -  0, 

(2c)           u(x,a)   -  0, 


-00  <  x  <  oo; 
-00  <  X  <  ooj 
-00  <  X  <  0, 


-  u  - 


2      2      2 

Here  k  (y)  =  o)  e(y)n/c  ,  and  c  is  the  velocity  of  light. 

We  will  not  specify  the  function  u  further  now.  However,  we  note  that  as 
a  consequence  of  the  symmetry  about  the  plane  y  =  a,  we  may  simplify  our  problem 
by  "bisection"  ,  i,e,,  we  can  replace  the  interval  0  <  y  <  2a  by  the  half  inter- 
val 0  <  y  <  a,  and  replace  condition  (2b)  by  the  condition 


(2c') 


au 
ay 


(x,a) 


0, 


0  <  X  <  00, 


The  relevant  situation  is  illustrated  in  Fig,  2, 
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Figure  2 


3»   A  priori  discussion  of  properties  of  the  solution 

In  this  section  we  shall  discuss  certain  a  priori  properties  of  the  solutions 
of  the  system  of  equations  (l)  and  (2).  Most  of  these  properties  follow  from  Sturm- 
Liouville  theory,  and  will  be  called  into  play  in  connection  with  the  developments 
of  Section  U, 

First  of  all,  we  form  the  product  solution 

u(x,y)  -  e^^  SZ(y) 
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for  Eq,    (1),   and  we  arrive   at  the  ordinary  differential  equation 
,2, 


(3) 


^  ^   [-  V^  *  k^Cy)]  q>   -0, 


2 
where  the  initial  condition  to  be  satisfied  by  0  =  0(-v  ,y)  is 

(U)  0(-v^,O)  -  0, 

and  where  the  function  is  normalized  by  a  condition  on  its  derivatives,  viz,, 

(5)  ?)'(-v^,C)  -  1. 

2 
The  function  ^(-v  ,y)  is  uniquely  defined  as  a  function  of  y  by  the  differ- 

2 
ential  equation   (3)  and  the  initial  conditions   (U),   (5),     Now,  writing  -v     »  A. 

for  convenience,  we  define  two  problems;   that  of  finding  eigenvalues  X  ■  ^i     such 

that 

(6)  ^{\^^,s)     "     0 

and  that  of  finding  eigenvalues  X  «  X  such  that 

(7)  ?l'(X^,a)  "  0. 

We  shall  refer  to  the  problem  of  finding  the  functions  <^{y.   ,y)  which  satisfy 
Eqs.  (3),  (U),  (5),  (6)  as  Problem  I,  and  that  of  finding  functions  (?(X  ,  y) 
satisfying  Eqs.  (3J»  (U),  (5),  (7)  as  Problem  II;  and  we  shall  denote  the 
solutions  ^{\i   ,y)  of  Problem  I  by  Y  (y),  and  the  solutions  ^(X  ,y)  of  Problem 
II  by  ^   (y).  Of  course  all  the  eigenvalues  in  question  are  real  and  simple. 

Notice  nov  that  the  eigenf unctions  T  (y),  ^  (y)  both  form  complete  sets  over 
the  interval  0  <  y  <  a.  The  function  u(x,y)  can  accordingly  be  expanded  in  the 
form 
(8a)  u(x,y)  -  J^     'n^^^^n^^^'       x  >  0; 

(8b)  u(x,y)  -  21  \^''^\^y^»  ^<°» 

where 
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(9a)*  vjx)   =  c^^^  exp(-  ^  x)     +     d^^^  exp(+  ^  x) , 


(2)  (2) 

(9b)  w^(x)  -  c^       exp(  VJr  x)       +     d^       exp(-  /iT  x) . 

The  expansions  (8a),   (8b)  are  obtained  by  applying  finite  integral  transforms  in- 
volving the  eigenfunctions  ^(y),  T  (y),  respectively,  to  Eq.   (i)   and  making  use 
of  the  appropriate  boundary  conditions. 

Next,  we  may  obtain  some  conclusions  about  the  eigenvalues  X  ,  u  .     Consider 
the  Stunn-Liouville  system 

(10)  ^"    +     fl  +  k^(y)1^    =    0, 

(11)  9  (>^,0)     -     '^(X,2a)     =     0 

over  the  interval  0  <  y  <  2aj  we   shall  call  this  Problem  III.     Notice  that   (10 ) 
governs  the  modes  in  the  full  guide   (for  x  >  O),  i.e.,  the  unbisected  guide.     Since 
K  (y)  >  0,  -we  see  that  there  may  be  negative  eigenvalues  X,  j  and  for  such  eigenvalues 
X  =  X^,  the  square  root  of  X  is  imaginary,  so  that  we  have   a  propagating  mode.     It 
follows,  then,  that  Problem  III  has  a  single  negative  eigenvalue,   say  X  ,   and  that 
all  other  eigenvalues  \    of  the  problem  are  positive. 

We  obtain  further  information  from  the   assumption  that  k(y)  is  even  about 
y  =  a.  Let  y  =  a+s,  where  -a  <  s  <  a;  then  9(X,a-s)   is  a  solution  of  our  problem 
if  p(X,a+s)   is  a  solution.     But  the  eigenvalues  are  assumed  to  be  simple,   so  that 
we  must  have  '^(X,a-s)   =  c  5f(X,a+s),  where  c  is    a  constant.   Setting  s  »  -s  in  this 
equation  we  deduce  that  ^(X,a+s)  =  c  ^(X,a+s),  or  c     «  1  and  c  =  t  1,  from  Which 
we  see  that  the  eigenfunctions  of  our  problem  are  either  even  or  odd. 

Now  arrange  the  eigenvalues  X.    in  increasing  order: 
\    <0<X^<\2<...      , 

\riiere  the  zero  is  interpolated  to  indicate  that  ^    is  the  only  negative  eigenvalue. 


For  r  =0,  we  write 


■„       -.i''exp(i/rrx).d^l'   _,..      „ 


v.  ■  cl"'   exp(i  V-V  *^  "^  %       exp(-i  ^-X    x) . 
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To  each  eigenvalue   A.     there  corresponds  an  eigenf unction  ^  (X.,y),  and  it  is  well 
known  (cf.    [3],  p.232)   that  ^  (X,,y)  vanishes     n     tiroes  within  the  range  0  <  y  <  2aj 
thus,  (?     (A.,y)  vanishes  2p  times  in  0  <  y  <  2a,     Notice  that  the  points  y  where 
9^-     vanishes  are  symmetrically  located  about  y  =  a,   in  virtue  of  the  even  (or  odd) 
symmetry  of  the  eigenf unction j  but  ^^     cannot  vanish  on  y  «  a.     If  it  did  have  a 
zero  at  y  ■  a,  this  would  have  to  be  a  double  zero   (otherwise  ^     could  not  vanish 
2p  times  in  the  interval),  and  a  double  zero  implies  the  vanishing  of  both  the 
function  and  its  derivative,  i.e.,  the  identical  vanishing  of  the  function.     By 
similar  reasoning,  it  follows  that  0p     ^   must,  on  the  contrsiry,  have  a  zero  at 
y  »  a,  to  preserve  the  required  odd  number  of  zeros  in  the  interval  0  <  y  <  2a» 
We  see  that  ^„     ,   is  an  odd  function  of  y,  while  ^_     is  an  even  function.     The 
latter  statement  implies  that  Pp   (0,a)  •=  0.     For  every  p,  then, 

(12a)  *$2p^^'°^     "     Kj>^''*^^     '     °5 

(12b)  ^2p+l^^'°^   "  ^2p+l^''»''^   "  °- 

Examining  Eqs.  (12a),  (l2b),  we  see  that  in  the  interval  0  <  y  <  a  the  even  eigen- 
functions  JL  are  of  the  type  ^  defined  above,  while  the  odd  eigenfxinctions  ^^ 
are  of  the  type  Yj  the  coiresponding  eigenvalues  are  Xp  ,  Xp  ,,  respectively.  On 
the  other  hand,  every  eigenfunction  of  type  ^   can  be  continued  by  even  reflection 
into  an  even  eigenfunction  of  the  double  interval,  while  every  eigenfunction  of 
type  T  can  be  continued  by  odd  reflection  into  an  odd  eigenfunction  of 
the  double  interval.  We  see  that  each  ^  is  of  type  ^p  ,  while  each  T 
is  of  type  ^2-p^i'   ^•®*' 

(13a)         gfp  -  ?2p,       Xp  -  >^p, 

_  ^  P  "  0,  1,  2,  ...; 

<"^>  ^p  ■  Vi'        ^P  ■  Vi 
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with  the  ordering 


\<°<^^o<\<^^< 


•  •  •      • 


This  is  as  we  should  expect:     the  n's  are  all  positive,  and  there  are  no  propagating 

modes  in  the  region  0<y<a,  x<0»     One  consequence  of  this  is  that  if  we  require 

a  solution  u(x,y)  which  is  finite   at  infinity,  we  must  have  d         ■  0   for  n  >  1,  and 

d         »  0  for  n  >  0  in  Eqs.   (9a)   and  (9b)   respectively, 
n  — 

A  further  a  priori  property  of  the  solution  as  it  is  so  far  prescribed  may  be 
obtained  by  a  mathematical  construction  related  to  energy  considerations.     We 
form 


/. 


3u     J 

U      •:?-       dS 

3n 


r  +  C 


over  the  contour  indicated  in  Fig.  3;  here    P  is  the  rectilinear  part  of  the  contour, 
ard  C     is  the  boundary  of  the  small  semicircle  of  radius  e  about  the  point  x  =  0, 
Then,  denoting  by  A     the  interior  of  the  contour,  we  have,   applying  successively 
Green's  formula  and  Parseval's  theorem, 


y  =  a 


y  =  0 


Figure  3 
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u§d..jj    |,u|W|u|ndS.(     u|    as 


o[_'o'  °Jl  °  °        ° 


(lU)  -  -i  X        |c^"M'  -   Id^-'M^U  J   2  X     Im  d^"^c^"^e 


rrrT-2iXox. 


,2 


-2X 


n^       .-        ,,    ,2  .,       ^^n^ 


-  7"  X    |c    r  M     e       "  '  -  ^"^l    Id    r  N     e 

a  a 

where  M^  -  |  |0^|^dy,  N^  -  j     ll^l^dy. 

^  o 

Comparing  the  real  and  imaginary  parts  of  (lU),  we  see  that 


lo<^'l'-M<^'i' 


d(i> 


and  since  the  reflection  coefficient  R  «  ,  t  ■  ,  we  have  found  that  R  ■  e  ^,  with 

% 
p  real.  This  result  was  to  be  expected.  Our  problem  is  to  find  the  phase  p.  We 


see  also  that  if  two  solutions  u-,  Up,  of  the  type  so  far  described  (i.e.,  expon- 
entially damped  except  for  an  incident  and  reflected  wave)  have  the  same  c   , 
then  they  have  the  same  d^   ,  For  consider  the  difference  v  -  u^  -  u„  of  these 
two  solutions;  v  is  a  wave  function,  and  we  have^ 

°o     °0,1   *^0,2'  %  %,1       %,2   • 

Since  c^^  -  0,  therefore  d^^  »  0. 
o      '  o 

Now  let  us  consider  two  solutions  u, ,  Up,  of  the  type  described  above,  both 

having  the  same  incident  mode  strength.  If  the  difference  v  of  these  solutions  is 
such  that 


Here  the  subscripts  1  and  2  are  used  to  indicate  the  coefficients  occuring  in 


iL  and  u  respectively, 


and  such  that 
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lim  I     V  -5—     ds 

e->0     ;         3" 


|\7  v|^   -  k^    |v|^l    dS 


exists,  where  A  is  the  region  0<y<a,  Xp<x<x^,  then  we  can  provB  that 

V  s  0»     Let.  V  be   any  wave  function  of  the  type  considered  here;    then  an  equation 

of  tte   same  nature  as   (lU)  holds  for  v»     Since  c^    '   »  0,  d         «  0  in  (Hi)  we  can 


o  '     o 

let  6  ->  0  and  obtain 


(15)  I    I  ^|Vy|2  -k^lvlHdS 


-  -  7ZV^   l^nl  V^L-2  V^  -l]    -  51  ^tjdj\exp[2  v^  x^] 


Here  c  ,  d  are  Fourier  coefficients  of  v,  and,  as  above, 
n'  n  >    *        » 

M„  -  j\^J  ^.  \     -  jiTj'dy. 

If  we  now  let  X- ,  x-  approach  zero,  the  double  integral  on  the  left  of  (i5) 
approaches  zero,  and  we  have 

00  —  00  r\ 

(16)  r"  X  ic  r  M    +  r  ti  Id  r  N    -  o. 

4—       n'n'       n  4—  ^n'n'       n 

Since  all  the  coefficients  in  (16)   are  positive,  it  follows  that  v  ■  0, 

We  may  interpose  here  the  remark  that  we  shall  find  a  solution    u     which  can 
be   shown  to  be  0^/p  sin   ((?/2)^  ,  where  p,  g  are  polar  coordinates  with 
pole  atx  =  0,y°a  and  where  ^  =  Oaty  =  a>.x<0,  and  0  <  {?  <  2n,     The  behavior 
of  7u  is  then  0{l/\/p),     This  behavior  of     u     and     Vu     can  be  shown  to  ensure  the 
fxilfillment  of  the  conditions  which  we  have  employed  in  sketching  our  uniqueness 

argument . ___^_^________ 

»  This  follows  from  the  previous  paragraph. 


-  n  - 


Continuing  our  review  of  general  properties  of  the  solution  of  the  problem 
we  have  posed,  we  notice  that  regardless  of  the  precise  nature  of  the  function 
k(y)  (provided  it  is  continuous),  we  have  for  Problem  III,  (cf.  [h]) 

2 


an  T    .    i^\    .  0(1) 


and  thus 


as  n  ->oo 


(18)  /T       =     i^^    +     0(i)  as  n  ^00   . 

RecalDing  that 


we  have 


(19a)  v^=     iigl^    .     0(i) 


as  p  — >0Oj 


(19b)  /IT  =     -^£ji^      *     0(i)  as  p  ^oo. 

pa  F 

In  what  follows  we  shall  also  need  the  asymptotic  expression  (cf.    [5],  p.  10) 
for  the  function  (?(X,y)   =  {2l(-v  ,y)  of  Problem  I.     The  result  in  question  is   that 


(20) 


;,(.v2,y)  ^  £i£iL^    [l.    0(i)]  , 


(21)  gf'(-v^,y)-cosh  vy     ll  +    0(i)l 

uniformly  in  y  (O  <  y  <  a) ,  as  |v|  ->oo. 

Eqs,  (20),  (21),  together  with  the  asymptotic  eigenvalue  formulas  (l9a), 
(l'^),  express  the  close  connection  between  general  solutions  of  a  Sturm-Liouville 
problem  and  the  trigonometric  functions.  In  addition  we  recall  (cf,  [3j )  the  well- 
known  result  that  0(X,y)  is  an  entire  function  of  X, 

Now  in  view  of  the  form  of  our  solution,  our  problem  reduces  to  that  of  deter- 
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mining  the  coefficients  c   ,  d^  ,   and  c^  '   ,  d^  '  '  for  n  >  1,  These  coeffi- 

oonn  — 

cients  are  determined  oy  an  infinite   set  of  equations  which  could  for  example  be 
deduced  from  the  condition  that  u  and  3u/3x  be  continuous  at  x  =  0,  by  comparing 
the   series  expansions  of  u  and  3u/3x  as  obtained  from  the  forms  which  are  valid 
for  X  >  0  and  x  <  0  respectively.     We   shall,  hovrever,    find  it  more  convenient  to 
follow  a  procedure  which  is,    as  we  shall  see,   closely  allied  to  the  function- 
theoretic  principle  underlying  the  Wiener-Hopf  procedure. 

We  can  present  our  solution,   for  all  x,   in  the  form  of  a  Fourier  integral, 
as  is  suggested  by  separation  of  variables,   and  as  can  be  deduced  by  the  applica- 
tion of  the  Fourier  transform  theorem  to  the  differential  equation.     This  repre- 

ivx 
sentation  is  based  on  the  improper  eigenfunctions  e         and  is   a  counterpart  to  the 

representation  in  infinite  series  which  we  deduced  from  the  finite  Fourier  trans- 
formations based  on  the  transverse  eigenfunctions  (jl  (y)   and  T  (y)  respectively. 
Our  procedure  is  to  set 
(22)  u(x,y)     -     ^     I   e^'^0(-v^y)   A(v)dv, 


where  C  is  a  countour  which  follows  the  real  v-axis  but  passes  under  the  real  poles 
of  A(v).     Experience  with  alternate  representations  of  solutions  of  the  wave 
equation  naturally  leads  to  the  supposition  that  the  discrete  expansions  which 
hold  for  X  >  0  and  x  <  0  are   simply  the  residue  expansion  of  the  integral  ex- 
pression for  u(x,y).     Since  exp(ivx)  is  exponentially  small  in  the  upper  half- 
plane  for  X  >  0   (and  in  the  lower  half-plane  for  x  <  0),  we  must  expect  the 
residue  expansion  for  x  >  0  to  be   associated  with  the  upper  half -plane,   and 
hence  we  expect  poles  of  A(v)   at  the  points  i  yx"  in  the  upper  half-plane, 
at  the  points  ±  >/-\     on  the  real  axis,  and  at  the  points  -i  •/yT  in  the  lower 
half -plane.     Then  if  A(v)^(-v  ,y)  is  of  algebraic  growth  at  infinity,  the  fact 
that  9f(-v   ,y)  is  entire  will  ensure  that  a  residue  expansion  is  possible,   and 

the  expansion  will  automatically  lead  to  series  representations  of  the  type  we 
*For  the  case  of  a  uniform  medium,   the   solution  of  the  infinite  set  of  equations  is 
actually  obtained  by  Hurd  I-  J ,  by  a  simplified  function  theoretic  method. 


-  13- 

we  have  discussed  above.     The  series  representations  will  ensure  that  the  necessary 
boundary  conditions  are   satisfied,   and  the  poles  on  the  real  axis  will  yield  the 
incident   and  reflected  propagating^  waves.     The  ratios  of  the  residues     will  provide 
the  reflection  coefficient.     Of  course,  to  carry  out  the   above  program  we  must  be 
sure  that  A(v)  has  no  singularities  other  than  the  desired  poles. 

The  fact  that  A(v)  must  possess  the  indicated  poles  can  be  directly  inferred 
from  the  seiles  expansions  if  we  apply  the  Fourier  transform.     We  have,  for  x  >  0: 

u(x,y)  o  c^^^expji  /IT  xj   +  d^^Wj^i  V^  ^  +  F  c^^^exp  [1  yTJ  xj(2(^(y) 


and  for  x  <  0: 


u(x,y)     -     y    c^^'  expj  v/IT  x    T^(y). 


Then  defining 


v(x,y)     »     u(x,y)   -  c^   ^exp  i  v/^  x     -  d^"*" ^exp  -i /TT"  x  ,       x  >  0, 
^^^^  -     u(x,y),  X  <  0, 

and  recalling  that  Y^(y)   =  (^{[i.^,y\9j^{y)  '  <?(\>y) ,  vre  have 

(2U)  I       V  e  dx  -  )_      +    2^    e 

<-oo  '"^n  '  n 

The  partial  fraction  expansion  (2U)  exhibits  the  imaginary  poles  explicitly  and 

shows  that  the  c^  '   and  c    are  associated  with  the  residues  of  A(v).  We  have 
n       n 

dealt  here  with  the  function  v,  instead  of  u,  since  the  Fourier  transform  of  the 
propagating  modes  does  not  exist  for  real  v.  From  this  very  fact  however,  we 
infer  that  A(v)  must  have  poles  on  the  real  axis  at  the  points  ±   \  , 

The  actual  construction  of  A(v)  is  quite  simple.  We  represent  it  as  a  fraction 
whose  numerator  is  an  entire  function  e(v)  and  whose  denominator  has  zeros  at  the 
■«•  at  the  ireaL  poles. 
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double  infinity  of  its  poles.     For  this  purpose  vre  define 

00 


(25) 


^^'^'-  -^  ?.(--^Hr^] 


p  /       L_    '  p 


(26)  JYJv)     =     ff    /l  +     -^  \    exp 

^  p=o  1,  i  VJT 


and  we  write 


(27)  A(v)  '  ^^"^ 


W^^p, 


1  + 


/IT" 


TTi(v)  T2(v) 


The  uniform  convergence  of  the  infinite  products  follows  from  the  asymptotic 
formulas  (I9a)  and  (I9b)  for  /JT*,  y^^  j  these  equations  demonstrate  the  con- 
vergence  of  the  two  series 


The  function  e(v)  must  be  an  entire  function,  and  we  shall  specify  it  later.  We 

2         2 
have  seen  the  analogy  between  the  fxinctions  9f(-v  ,y),  (2'(-v  ,y)  and  the  sine  and 

cosine  functions  respectively:  we  now  call  attention  to  the  anaD.ogy  between  lTi('^)> 

7r«(v),  and  the  gamma  function.  This  analogy  is  suggested  by  the  forms  of  the 

products  (25),  (26)  and  also  by  the  asymptotic  character  of  their  zeros}  it  is 

made  explicit  if  we  consider  the  f\mctions 

(28)      F^  .  TTi(v)  Tri(-v)  -  F^(v,a)  -  A  +  ^  \  Tf  (l  *  I"  ) 


00  /    2  \ 

Fg  -  TT2(v)  Tr2(-v)  -  F2(v,a)  -   fT^  ^1  +  ^-  j  • 


The  dependence  of  the  functions  F, ,  Fjj  on  the  quantity  a  occurs  via  the  eigen- 
values X  ,  u,  ,  which  are  the  roots  of  the  equations  g*{\   .a)  »  0,  p(|a  .a)  -  0, 
n'  '^n'  n  n 


respectively. 
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Now  consider  the  asymptotic  formulas  for  the  entire  functions  Qf(-v  ,a), 
(^'(-v  ,a).  From  these  formulas  we  see  that  the  functions  are  of  order  1;  and 
the  zeros  are  in  arithmetic  progression,  so  that  their  exponent  of  converf^ence 
is  also  1.  Hence  the  genus  of  (fj   and  0'  is  1,  and  we  have  from  the  Hadamard 
factorization  theorem  (cf.  [6],  p.  250 ) 


(30)    iZ(-v^a)  -  9f(0,a)  TT  (  1  " 


exp 


or 


(31)    5>(-v^a)  -  9((0,a)  Tr2(^)  TT2(-^)  '  (^^^'^^  F2(v,a). 


Next,  consider  the  asymptotic  formula 


r4(   V.2   X    sinh  V  y 
7(-v  ,y)  -~ — ^ 


For  the  latter  function  we  have  the  identity 

sinh  n  v  n 


^  r(l  +  iv)r(l  -  iv) 

Tt-  2 

Comparing  with  (31)  we  see  that  the  relation  of  the  function  TTo^'^)   "^°  9'(-v  ,a)  is 
analogous  to  the  relation  of  the  function  l/  fCl  ■•■  1^)  '^^  'the  trigonometric 
functions.     In  addition  to  the  product  representation  for  ^^'j  we  have  incidentally 


(32) 

or* 
(33) 


9i'(V,a) 


0'(o,a)/'i*f  ^fr/i 


-^ 


i  fT] 
I'  p  / 


exp 


V 


i  sTT 
^  P 


00  /  v^ 


I  exp 


-V 


i  JTT 

^  P 


{Z|.(-v  ,a)   -  j?'(0,a)  Tri(v)  ITi(-v)   -  F^(v,a). 


It  is  useful  to  investigate  the  growth  of  the  products  Jf,  (v),  "[Tc(''')»     ^® 

quite  naturally  compare  them  with  the  gamma  function.     Since  we  have 

It  is  interesting  to  note  that  we  have  the  identity  (/(-v  ,2a)^(-v  ,a)9f'(-v  ,a) 
gf(0,2a)/v((0,a)(;i'(0,a),  which  is  a  generalization  of  the  equation 

y^Zil    /3in>^?:7y  ^^^7       „  2 


sin 


v(?w 
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VxT  -(p  *  |)  7   *  o(^)  ,  (p  -  0,  1,  ...), 


p        •'        2     a  -    p 


y;P  ^(p  4  1)  ^     +     0(|), 


(p  =  0,  1,    ...), 


we  consider  the  function 


00    / 


(3U,      ^(.iJV).^A,_^v      N    ^/ 


av 


r(|) 


1  \         ni(p  *  2^  J 


Tii(p  +  ^) 


ttH'-^H 


where  Y  is  the  logarithitdc  derivative  of  the  gamma  function. 
Then  in  accordance  with  (3U), 


(35) 


TTi(-v) 


1  + 


^J  I   U-^av/i(p+  ^)n 


exp<-  J 


^fl  1 


J\K    (p^?> 


r) 


In  deriving  (35)  a  rearrangement  has  been  undertaken  which  is  permissible  in  virtte 
of  the  separate  convergence  of  the  various  infjjiite  products;  the  procedure  here 
and  in  what  follows  is  the  same  as  that  employed  in  a  similar  situation  i-  -■,  We 
rewrite  our  expression  (35)  in  the  form 

Tri(-v) 


(36) 


V^^) 


1  + 


v^ 


p^  P(v)  exp(a-j^v) 


where  a-  is  a  constant  given  by  the  series 


"^'^?{^-^?)- 


Because  of  the  asymptotic  form  of  v^  >  this  series  converges;  the  constant 


Pi   -  TT 
■^        1 


00   (p  +  j-)- 


a 


^  p 


is  similarly  seen  to  converge.  The  function  P(v),  defined  by 
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P(v) 


00  i(  V^  -  T-.(p  +  i)/a)  /v 

IT  ^1^   — ^x ^-~-    ^ 

i  (p  +  K-)n 

i  £ —  .  i  +  1 

\^  a  V  j 


can  easily  be  shown  to  converge  unifonnly  as    |v|  ->oo  in  the  upper  half  of  the 
v-planej  the  method  is  again  the   saine   as  that  used  in    [2],  and  depends  on  the 
asymptotic  form  of  the  /T^*     Then,  letting    |v|  ->ao  in  the  upper  half-plane,  we 
get  P(v)  -^  1,   and 


(37)       lim 


iav       3 


lim 


^''^;!;i;^^^'  -[¥  4)]  ■  (-  ^)  ^.-(  V)- 


In  a  similar  way  we  introduce  the  function 

00 


and  we  find  that  as  v  — >oo  in  the  upper  half-plane 


(39) 


where 


lim 


P3exp(a3V)  Ul.-iL-^ 


".-fr^ 


\V^  (P 


+l)Tt/a  ) 


oo 

Pc  -  IT 


(p  •fl)Ti/a 


We  have  defined 


(27) 


A(v) 


e(v) 


1  + 


'      o 


Tri(v)Tr2(v) 
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Hence,  our  results  on  9^(-v  ,a)  permit  us  to  write  also 


A(v) 


e(v)  TToC-^) 


i*-^^lTi(v)  TTjCv)  ITjC-v) 


or 


But  we  have 


e(v)  TToC-v) 


A(v) 


^(0,a) 


(-^)Tr,(v)   <«--' 


e(v)      'H'l^-'^^    i?'(0,a) 
\   TToCv)   9'(-v^a) 


1  + 


'   o 


TTiC-v)  TTiC-v) 


lav^s 


^2  (")    1  «L<-  ^' 


QgC-  iav/n)  1  Q^(-  iav/n) 


and  hence  we  find  that  as  v  ->oo  in  the  upper  half -plane. 


TTi(-v) 

lim  ■  lim 

TTs  (V) 


ivaTT^ 


r(i)  r(-^-|)  exp(i^T(4)) 


(UO) 


=     0(vV2)     ^^(i^) 


(ou-ao)n 


Y(2)  -  T(|)   *       ^     ^ 


ai 


The  last  estimate  is  obtained  by  applying  Stirling's  formula  to  the  quotient  of 
gamma  functions;  we  shall  have  occasion  to  use  this  estimate  later.     For  convenience 
we  now  introduce  the  notation 

(Ul)  T(2)  -  Y(|)  +  (5^-^)    n     -     q. 
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2 

Now  we  shall  examine  the  growth  of  the  function  A(v)«  (^{~v   ,y) ,  Using  the 

2 
asymptotic  formula  for  (ji(-v  ,y)   we  find,  in  the  upper  half -plane, 

1  3 

(U2)  A(v)9(-v^y)  .  e(v)  0(v"  ^)exp(i^)  ^^^^   "   °(^"  ^)e(v)exp(i^) . 

In  the  lower  half-plane,  we  employ  the  expression 


(^3)  ,(v)     .     e(v)  Iti^    J12^ 

J[^   (v)     9f(-v%a) 


to  obtain  the  same  result.     To  ensure  algebraic  growth  at  infinity     we  choose 

(Ui)  e(v)     »     exp(^-    ^^  . 

Although  it  is  not  our  purpose  to  discuss  the  question  of  uniqueness  per  se  in 
the  present  work,  a  few  remarks  are  called  for  here  in  connection  with  the  choice 
of  e(v).     The  present  choice  of  e(v)  is  not  uniquely  determined.     We  might,  in- 
stead, have  taken 

e(v)     =     V  exp  ^ —^  \   . 

This  would  correspond  to  taking  the  x-derivative  of  the  solution.  Now  our  solu- 
tion in  the  present  form  can  be  shown  to  be  0(x  '  )  on  the  line  y  «  a,  as  x  ->  0. 

-1/2 
Multiplication  of  e(v)  by  v  would  yield  a  solution  which  would  be  0(x  '  ),  Such 

a  behavior  would  be  inadmissible  on  physical  grounds.  On  the  other  hand  our  solu- 
tion  is,  except  for  a  constant  factor,  the  only  solution  which  is  0(x  '    )   ,     The 
constant  factor  referred  to  here  arises  becaase  we  have  implicitly  prescribed  the 
complex  amplitude  of  the  incident  mode  in  our  solution.  It  is  in  fact  the  re- 
sidue of  A(v)  at  V  -  +  V-X  ■  .  To  obtain  a  different  constant  amplitude  we  must 
multiply  the  solution  by  a  constant  factor.  However,  the  reflection  coefficient 

*  The  proof  depends  on  the  fact  that  the  reqiiirements  indicated  in  our  earlier 
discussion  of  uniqueness  are  fulfilled  for  such  solutions. 
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is  independent  of  this  factorj  it  is  given  by 

(U5)  R  . 


!fJhtfl?(-i)-($)-(-^ 


qa 


00 


Since  /-X  ~  and  q  are  real,  it  is  evident  that  |R|  is  unity. 
Let  us  now  recapitulate  our  solution.  It  is 

(U6) 


u  -  ^   /  ®^'"'  A(v)0(.v^y)  dv. 


whei^ 


(U7)      <?"(-v^,y)  +  (-v^  +k^(y))  0  -  0;   9((-v^,0)  -  0}   (2)'(-v^0)  -  1, 


and  where 


|exp(-ivaq/n)J 
(U8)    A(v) 


in  this  last  equation  \  ,  [i     are  the  respective  sets  of  roots  of  the  equations 
9f'(\  ,a)  •  0  and  (jl{\i.   ,a)  =  0.  All  these  roots  are  positive  except  \^,  which  is 
negative.  The  contour  C  in  Eq.  (U6)  runs  along  the  real  axis  except  for  indenta- 
tions below  the  points  v  «  ±  /-\  ".  The  constant  q  of  Eq.  (U8)  is  given  by 


o 


q  «  T(2)  -  T(|)  + 


^|f  (v^^P^y  "^y^"  (P^i)j)| 


The  calculation  of  A(v)  is  considerably  simplified  if  we  express  the  exponential 
in  the  numerator  as  a  product  of  exponentials  in  accordance  with  the  definition 
of  q  as  a  sum,  and  then  divide  each  convergence  factor  in  the  denominator  by  the 
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corresponding  factor  of  exp(-ivaq/n) ,     Then  we   get 

exp[-^(T(2)-T(|))] 
A(v)    -      ^ — — , 

E  *  ^]  F^^"'/^  ^^^"^  r^^^P*  i^i)]  1^(1*^/^  /ir)exp(iv/(p+i)|)J 

We  do  not  propose  to  discuss  the  verification  of  the  solution  in  this  article, 
since  our  main  purpose  is  to  discuss  the  relation  of  the  solution  to  Sturm-Li ouville 
theory.     Suffice  it  to  say  that  the  verification  can  be  carried  out  by  standard 
methods  of  verifying  contour  integral  solutions  of  partial  differential  equations, 
once  we  are  assured,   as  we  already  are,  that  the  function  A(v)  ^'(-v  ,a)  is  of 
algebraic  growth.     The  method  is  illustrated  in  detail  in    [2],  and  depends  on 
deformation  of  the  contour  of  integration  so  as  to  allow  differentiation  under 
the  integral  sign. 

U.       Concluding  Remarks 

It  should  be  noted  that  the  problem  treated  here  is  the  simplest  example  of 
a  class  of  problems  amenable  to  the  present  treatment.     In  the   general  case  we 
wDvild  have  an  asymmetrically  placed  half-plane,   asymmetrically  varying  (or  even 
discontinuously  varying)  t  and  ti,   mixed  boundary  conditions  on  all  the  conductors, 
and  an  jarb^traiy  number  of  propagated  modes.     We  would  then  have  to  deal  with  the 
Sturm-Liouville  theory  of  three  regions  of  space,  but  nevertheless     the  theory 
developed  here  can  be  carried  over  on  the  lines  indicated  above.     The  lengthier 
Wiener-Hopf  integral -equation  procedure  could  be  employed  in  all  such  problems. 
The  use  of  that  method  ultimately  reduces  to  the  factorization  of  the  Sturm-Liouville 
functions  and  their  derivatives  as  functions  of  the  transform  variable.     In  this 
context  the  present  work  indicates  a  general  and  convenient  solution  of  that  factor- 
ization problem.     In  the  course  of  the  work  we  have  supplied  a  factorization  of 

2  2 

the  functions  9r(-v  ,a)  and  ^»(-v  ,a).   [See  also  the  footnote  on  page  1$Q  In 

special  cases,  the  results  obtained  by  our  method  of  course  reduce  to  the  results 
previously  obtained  by  other  methods  by  various  authors. 
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12.  Abstract: 

A  simple  solution  of  a  general  problem  involving  a  bifurcated  waveguide  is 
presented.  The  purpose  of  the  work  is  to  explain  a  new  and  simple  method 
of  solving  such  problems  and  to  exhibit  an  organic  connection  between  Sturm- 
Liouville  theory  and  the  theory  of  two-part  boundary- value  problems. 


